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We investigate /(i?)-gravity models performing the ADM-sIicing of standard General Relativity. 
We extract the static, spherically-symmetric vacuum solutions in the general case, which correspond 
to either Schwarzschild de-Sitter or Schwarzschild anti-de-Sitter ones. Additionally, we study the 
cosmological evolution of a homogeneous and isotropic universe, which is governed by an algebraic 
and not a differential equation. We show that the universe admits solutions corresponding to 
acceleration at late cosmological epochs, without the need of fine-tuning the model-parameters or 
the initial conditions. 
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I. INTRODUCTION 

Recently, there has been a tremendous thrust in re- 
search activities in cosmology, due to the wealth of data 
from various experiments that are already available and 
more that are anticipated in the near future. Amongst 
others, the discovery of the late time acceleration of the 
universe has been particularly fascinating for cosmolo- 
gists, particle physicists and string theorists alike. In par- 
ticular, the significant improvement of the data-statistics 
(from 50 to more than 300 Supernova(Snla)) [l[ has made 
the aforementioned result indisputable. Furthermore, 
complementary probes like Baryon Acoustic Oscillation 
01 and Cosmic Microwave Background Radiation mea- 
surements, [1] offer additional evidences. 

Although the simplest candidate of the acceleration 
mechanism is the cosmological constant, one can con- 
struct various "held" models in order to incorporate this 
mysterious "dark energy" In particular, one can use a 
canonical scalar field (quintessence) Q , a scalar held with 
non-standard kinetic term (k-essence) [1] or with a neg- 
ative sign of the kinetic term (phantom) Q , the combi- 
nation of quintessence and phantom in a unified model 
named quintom [3], scalar helds non- minimally coupled 
to gravity Q or simple barotropic fluids with specific 
pressure-form such as Chaplygin gas [9| (for nice reviews 
on dark energy models, see [lO|). 

Instead of using field dark energy constructions, an 
alternative approach is to modify gravity itself. The 
Dvali-Gabbadadge-Poratti (DGP) the Cardassian 

[T^ and the Shtanov-Sahni [3] models follow this di- 
rection, lying in particular in the higher-dimensional sce- 
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nario sub-class. However, instead of using extra dimen- 
sions, one can insert higher-order curvature invariants in 
the usual Einstein-Hilbert action, with the simple con- 
sideration the so-called /(i?)-gravity models that 
is the addition of Ricci-scalar functions. Such mod- 
els wish to alleviate the non-renormalizability of grav- 
ity, and acquire theoretical justification from low-energy 
string theory. Alternatively, since higher-order terms 
can be related to non-minimally coupled scalar degrees 
of freedom, these models are equivalent to scalar-tensor 
constructions. Concerning their cosmological implica- 
tions, they can describe both inflation as well as the 
late-time acceleration. However, most of /(i?)-gravity 
models do not manage to pass the observational and 
theoretical tests (solar system, neutron stars and bi- 
nary pulsar constraints), giving also rise to an unusual 
matter dominated epoch and leading to signiflcant flne- 
tunings [Tsj . Moreover, even if improved versions (like 
the "Chameleon mechanism" [l6i] ) manage to satisfy the 
above constraints, one can still face problems in the 
strong gravity regime, due to the curvature singularity 
that is inevitable at the nonlinear level p7j . 

Although higher time-derivatives can be beneficial in 
making gravity renormalizable, they also lead to ghosts. 
However, the recently developed Hofava gravity wishes 
to act as a power-counting renormalizable, Ultra- Violet 
(UV) complete theory of gravity without possess- 
ing the full diffeomorphism invariance of General Rela- 
tivity but only a subset that is manifest in the Arnowitt, 
Deser and Misner (ADM) slicing. There has been a large 
amount of effort in examining and extending its prop- 
erties, as well as exploring its cosmological implications 

Motivated by these, in the present work we are in- 
terested in investigating f{R) gravity models performing 
the ADM slicing of standard General Relativity, that is 
its (3-|-l)-decomposition based on the Hamiltonian for- 
mulation [20| . In particular, we wish to extract the 
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static, spherically-symmetric vacuum solution for general 
/(i?)-models under ADM decomposition, and study the 
homogeneous and isotropic cosmological solutions which 
present late-time acceleration. The paper is organized 
as follows: In section |lT] we present a brief introduction 
on ADM (3-|-l)-decomposition and we write the gravi- 
tational action for /(i?)-gravity models. In section Hill 
we derive the static, spherically-symmetric vacuum so- 
lutions for general f{R) models under ADM slicing. In 
section HV] we apply this approach to cosmological frame- 
works, investigating universe evolutions that experience 
late-time acceleration. Finally, section |V] is devoted to 
the summary of the obtained results. 



III. STATIC SPHERICALLY SYMMETRIC 
SOLUTIONS 

We are looking for the static, spherically-symmetric 
vacuum solutions of the aforementioned general f{R) 
gravity models, under ADM decomposition. In this case 
the metric ([ij writes: 



fir 



(7) 



that is Kij = and ^ _^[r/'(r) + /(r)-l]. Setting 
f{r) - I = X{r), we obtain ^/'(r) + f - I ^ rX'{r) + 
X{r) = [X{r)r]'. Thus, defining rL{r) as 



II. ADM (3+l)-DECOMPOSITION AND f{R) 
GRAVITY MODELS 

We start by considerin g th e ADM decomposition of the 
four dimensional metric 12011: 



ds^ 



-N^df + gij{dx^ + N'dt){dx^ + N^dt). (1) 



Here N and are the lapse function and shift vectors 
respectively, while gij is the induced metric on the 3- 
dimensional space-like hypersurface for a fixed time. In- 
dices of all projected tensors can be lowered (raised) by 

For the space-like hypersurface with fixed time, extrin- 
sic curvature is defined as 
1 



K, 



2N 



(2) 



where dot represents derivative with respect to time and 
i,j = 1,2,3. Under such decomposition. 



(3) 



where i?*^^^ and R'^'^^ are the four and three-dimensional 
Ricci scalars respectively. Hence, the Einstein-Hilbert 
action can now be written as 
1 



S = 



WnG 



dtd^xNy^{R^^^ 



where g'^ is the determinant of the three-dimensional 
space-like hypersurface, and Sm accounts for the matter 
content of the universe. Generalizing the above action, 
one can write the action for an /(i?)-gravity model as 

j j^dtd^xN^f{R^^^+K,,K'^-K^)+S„,. 

(5) 

The action is in general function of N, A^* and gij and 
hence varying the action with respect to these quantities 
leads to the equations of motion. In our subsequent cal- 
culations we assume /(i?*'')) = R'^^^ + thus the 
above action becomes. 
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IGttG 



dtd^xNy^ + K,jK'^ - K^) 

- S,n- (6) 



R^'^^--[X{r)r]'^n{r), 
the action ([6]) after angular integration reads 



(8) 



S ^ 



1 



IGttG 



dtd^x 



VW) 



[H{r)+F{H{r))]. (9) 



Varying ([9]) with respect to N{r) and setting SS/SN{r) 
we obtain 



nir)+F{nir)) = 0. 



(10) 



This is an algebraic equation for 7i(r) depending on the 
functional form of F. One can always solve this equation 
to get 



7Y(r) ~ constant. 



(11) 



which depends upon the functional form of F{TL{r)). De- 
noting the above constant by the parameter /3, we can 
obtain 



H{r) = -^[X{T)r]' = 13 



X{r) 



A 
r 



(12) 



where A is the integration constant, set from now on to 
— 2M with M a new constant. Thus, 



r 

Variation of ^ with respect to f{r) leads to 
d / dL\ dL 



dr\df'J df 



0, 



where 



N{r)r 



.2 r 



V7 



nifJ',r)+F{nifJ',r) 



(13) 



(14) 



(15) 



A straightforward calculation gives 
dL 1 N{r)r^ 



df 



n + F{n) 



2 /3/2 

N{ry I an dF{n) 
~7T 



dR^ 



(16) 



which under (flOl) leads to 



by Starobinsky in [2l|: 



dL 
df 



2N{r) 



dF{H) 



(17) 



Similarly, we acquire 



N{ry p m dF{n) 

2N{r)r 



dF{n)- 



V7 



1 



(18) 



In conclusion, inserting (|16p . p8)) into p4|) . and using 
that 1 + ^f^j^J^)^ = constant (arising from the solution 
for f{r) in i?'^^)) finally obtain 



d_ 

dr 



N{r) 







N\r) = f{r). 



(19) 



This result, together with equation pS]) . shows that the 
static, spherically-symmetric vacuum solution for a gen- 
eral f{R) model under ADM decomposition is either a 
Schwarzschild de-Sitter or Schwarzschild anti-de-Sitter 
one, depending upon the choice of /3. 



fiR) = Ai?o 



- 1 



(21) 



where A, Rq and n are the model parameters, and from 
now on, R denotes the four dimensional Ricci scalar. The 
advantage of choice (PTjl is that it accepts a theoretical 
justification, but one could also use a different ansatz at 
will. Using this choice for f{R), and under the metric 
(|20)) . the action ([6]) becomes 



S = 



167rG 



- 1 



} + (22) 



where = - is the Hubble parameter. Finally, as usual, 
the energy density and pressure for the matter field are 
respectively defined as 



IV. HOMOGENEOUS AND ISOTROPIC 
COSMOLOGICAL EVOLUTION 

Here we assume that the background is homogeneous 
and isotropic and the spatial 3-hypersurface is flat (fc = 
0): 

ds^ = -N^{t)dt^ + a^{t){dx^ + dy^ + dz"^), (20) 

where a{t) is the scale factor and N{t) is the lapse func- 
tion. In order to present an example of explicit cosmo- 
logical solutions, we consider the form for f{R) proposed 



J 



1 5S„ 



9ijPm 




(23) 
(24) 



Variation of action ((22|) with respect to N, and the 
subsequent fixing = 1 (as it is usual in cosmologi- 
cal applications of the "foliation preserving" framework), 
leads to the "effective" Friedmann equation 



XRo r 



6 



24XnH^ 

Rn 



7?2 



-Pn 



(25) 



r 



According to the usual approach, one could also vary the 
action ((22|) with respect to the second variable gij. This 
procedure, although straightforward, leads to a compli- 
cated result, which forbids its physical use. Alternatively, 
we prefer to use the matter conservation equation: 



3H{prn +Prn) = 0, 



(26) 



the matter to be dust [pm = 0), which is definitely jus- 
tified for late-time cosmological behavior. In this case, 
(I^S)) leads to the usual evolution pm = Pmoo,~^, with pmo 
the matter energy density at present, where the scale fac- 
tor is fixed to 1. Therefore, inserting this formula into 
(l25l) we result to 



assuming for simplicity, and without loss of generality. 
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Hi 



\a r 
"6 



Rl 



2A\n / H 



r 



(27) 



where we have defined a = with iJo the current H- 
value and r^mo the matter density-parameter at present. 
Thus, the parameter a accounts for the modification of 
gravity. Finally, note that imposing (P?]) at present, that 
is taking a — 1 and H = Hq, allows for the elimination 
of the parameter A in favor of a, n and i^mo, namely: 



A = 



1 



+ ^(l + f) 



-(n+1) ■ 



(28) 



Equation (j27p is the modified Friedmann equation of 
the modified-gravity model at hand, and contains all 
the cosmological information of the system. It presents 
the significant advantage that it does not contain any 
higher-order time-derivatives that appear in the usual 
approach of /(i?)-gravity models. Instead, and due to 
ADM decomposition, it is just an algebraic equation for 
the Hubble-parameter H(a), although of not simple form. 
Therefore, one can examine its solutions for various val- 
ues of n, which in turn can be used to construct all the 
observable quantities like deceleration parameter, lumi- 
nosity distance, angular diameter distance etc. 

In particular, knowledge of H{a) allows for a straight- 
forward calculation oidH{a)/da = H'{a), while for every 
quantity Q we obtain Q = Q' {a)aH{a). Therefore, for 
the deceleration parameter q = —a/[aH{a)]'^ we acquire: 



H{a) 



H'{a). 



(29) 



As usual, q <Q corresponds to a > that is to an acceler- 
ating universe, while g > corresponds to a decelerating 
one. Finally, q < —\ corresponds to > 0, which is the 
case of a super-accelerating universe p^ . 

As it is usual for modified gravity models, the role of 
matter is crucial for the determination of the cosmolog- 
ical behavior. For instance, in the complete absence of 
matter, that is setting f2mO to zero, (1271) implies immedi- 
ately that H is independent of a, that is H' {a) = = H 
and q = —1 which is just what is expected in this case. 
In the following we focus on the realistic case where 

^niQ ~ 0.3. 

As can be seen from (|27l). even for n — 1 the equa- 
tion for H{a) is of high order, giving rise to many solu- 
tion branches, the number of which increases fast with 
increasing n. Some of these solution-branches lead to 
imaginary H{a) and thus are not physical. Addition- 
ally, one gets solution-branches that lead to divergencies 
in finite scale factors in the past, which must also be 
omitted. We are interested in those branches that have 
H{a) > at all a and q{a) < at large a, that is cor- 
responding to an expanding universe which accelerates 



at late cosmological epochs. One significant advantage 
of the aforementioned procedure is that since we do not 
solve any differential equation, we do not have to face 
the discussion of the initial-condition determination. We 
only have to fix the values of Hq, Qmo at present, and 
then the system is fully determined for particular values 
of n and a. 

Let us first investigate the n = 1 solution subclass. 
In fig. [T] we depict the deceleration parameter q{a) as it 
is given by (|29p for the numerically obtained solution of 
the algebraic equation (|?7|) for n = 1 with ^Imo = 0.3. 
The curves correspond to four values of the parameter a. 
As we observe, in all cases we do obtain acceleration at 




FIG. 1: The deceleration parameter q{a) for an expanding 
universe for n — 1, with Slmo = 0.3 and Ho = 1, ao = 1. 
The curves correspond to a = 0.1 (black-solid) , a = 1 (red- 
dashed), a — 5 (green- dotted) and q = 10 (blue-dashed- 
dotted). The horizontal line marks the q = bound. 



late times, with the transition to the accelerated phase 
realized at earlier times for larger a-values, that is for 
more significant gravity-modification. It is interesting to 
notice that for large values of a, the behavior of q{a) can 
be non-monotonic (as can be seen in the a = 10 curve), 
and going to even larger values (a = 20) it brings de- 
celeration at very late times. Since this scenario is not 
favored by observations we do not show it explicitly, but 
it is characteristic of the rich phenomenology and cos- 
mological possibilities that our model presents. Finally, 
we mention that in all cases q{a) is larger than —1, that 
is the imposed /(i?)-ansatz does not seem to be able to 
lead to a super-accelerated universe. 



In fig. [2] we depict the g-behavior for the n = 2 solution 
subclass. In this case, the curves for a — 0.1 and a = 1 
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FIG. 2: The deceleration parameter q{a) for an expanding 
universe for n = 2, with fimo = 0.3 and Ho — 1, ao — 1. 
The curves correspond to a = 0.1 (black-soUd) , a = 1 (red- 
dashed), a — 5 (green-dotted) and Q = 10 (blue-dashed- 
dotted). The horizontal line marks the q = bound. 



have a small difference (not observed in the scale of the 
figure), and one needs to go to larger a in order to see 
a different cosmological evolution (which is achieved fast 
for a becoming larger that 1). As we see, we again obtain 
acceleration at late cosmological epochs. Note that for 
a = 5, the transition to the acceleration phase is realized 
earlier than for a — 10, as a result of the highly non- 
linear behavior of equation (p7| . 

In fig. [3] we depict the g-behavior for the n = 3 solution 
subclass. Similarly to the previous cases, we do obtain 
late-time acceleration, with q{a) for large a being non- 
monotonic. Finally, we mention that qualitatively similar 
results arise for larger values of n too, but for simplicity 
we do not present them explicitly. 

The aforementioned solutions correspond to H{a) > 
at all a's, that is to an expanding universe. However, for 
completeness we mention that the present model allows 
also for solutions that describe a contracting universe. 
Indeed, since (|27|) is a even equation for H{a), we de- 
duce that for every 77(a)-solution, ~H{a) is also a solu- 
tion. Therefore, while the investigated solutions of fig- 
ures[T]to[2]possess H{a) > at all scale factors, the same 
g(a)-behavior arise from the corresponding branches with 
H(a) < at all a, that is a contracting universe. 



V. CONCLUSIONS 

In the present work we have studied /(i?)-gravity mod- 
els performing the ADM slicing of standard General Rel- 
ativity, that is its (3-|-l)-decomposition based on the 
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FIG. 3: The deceleration parameter q{a) for an expanding 
universe for n = 3, with SlmO = 0.3 and Ho = 1, ao = 1. 
The curves correspond to a = 0.1 (black- s oltd) , a = 1 (red- 
dashed), Q = 5 (green- dotted) and a = 10 (blue-dashed- 
dotted). The horizontal line marks the g = bound. 



Hamiltonian formulation. This approach allows for an 
easier treatment of modified-gravity systems and for the 
extraction of their general theoretical and cosmological 
implications. 

As a first application we derived the static, spherically- 
symmetric vacuum solutions for general /(i?)-ansatzes. 
As we saw, they correspond to either Schwarzschild de- 
Sitter or Schwarzschild anti-de-Sitter ones, depending 
upon the choice of a particular parameter (3. 

Concerning applications in a cosmological framework, 
we investigated the evolution of a homogeneous and 
isotropic flat universe. Imposing as a specific example 
a particular ansatz for f{R), we showed that the Hubble 
parameter is given by an algebraic equation in terms of 
the scale factor, with a new parameter that determines 
the modification of gravity. This fact is an advantage 
since one does not need to discuss the initial condition, 
since all the information is included in the value of the 
matter density parameter at present. The system accepts 
many solution branches, the physical sub-class of which 
corresponds to either expanding or contracting universes. 
Furthermore, one can easily acquire solutions that cor- 
respond to acceleration at late cosmological epochs, in 
agrement with observations, and this is achieved with- 
out the need of fine tuning the model-parameters or the 
initial conditions. The model at hand presents rich cos- 
mological behavior, and moreover one can calculate ad- 
ditional observables such as luminosity distance and an- 
gular diameter distance. 

ADM decomposition, which is crucial in General 
Relativity and its possible quantum-gravitational ex- 
tensions such as Hofava theory, proves very convenient 
in studying modified-gravity constructions in general. 
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One could apply this formalism in order to extract 
the black-hole properties, to study the gravitational 
wave spectrum, to examine the cosmological bounce, 
etc. Finally, one could use observational data in order 
to constraint the possible gravity modifications, both 
qualitatively and quantitatively. Such a research is left 
for a future investigation. 
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